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1 Introduction 

The rigid rotator is a well-known problem in classical mechanics 1 as also in quan- 
tum mechanics. The quantum theory of the band spectra of diatomic molecules began 
with the quantisation of angular momentum of the rigid rotator. Much of spectroscopic 
calculation relies on rigid rotator models. It is well-known that there exists a fundamen- 
tal difference between atomic spectra and molecular spectra. While for atoms the line 
separation in a series decreases rapidly, for molecular spectra in the infrared, it is ap- 
proximately a constant. Actually, for linear diatomic molecules, two additional modes of 
motion, namely rotation and vibration, not present in atoms, are possible. In literature, 
the energy spectrum of such a system is investigated non-relativistically. Presently, we 
propose to view the model in the relativistic framework. 

After the classic application of relativistic wave equations to the atomic spectrum of 
hydrogen atom, there were several other applications of the relativistic wave equations 
to well-known problems of non-relativistic quantum mechanics. In the simplest from, 
one can cite the work of Sauter 2 , who treated the problem of a particle in an uniform 
field relativistically. Srinivasa Rao et.al. 3 , have dealt with the relativistic rigid rotator 
in an elegant way, treating the rotor as a dumb bell with equal masses. Following the 
prescription of Rao, we extend the treatment to a linear rotator with unequal masses. 
The foregoing discussion is applicable to any heteronuclear system, special case of which 
is the homonuclear species. Such studies, apart from being pedagogical in nature are 
potentially exciting and shed light on the spectra of diatomic molecules. 

2 Klein- Gordon Rotor 



In the simplest model of a rotator, we consider two atoms of masses mi and m 2 
separated by distance r. 

The Klein-Gordon Hamiltonian for a system of two non-interacting particles may be 

written as , , 

H = cy mlc 2 + pi + cy mlc 2 + P2 (1) 

where p\ and p 2 refer to the momenta of the two particles. 

Introducing the centre of mass co-ordinate R and the relative co-ordinate r, the 
two body problem may be reduced to an equivalent one-body problem. The relative 
co-ordinate, r : (x, y, z) is expressed as 

r = T\ — r 2 (2) 



and the centre of mass co-ordinate R : (X, Y, Z) is given by 



R= m ^ + (3) 
mi + m 2 



It is straightforward to verify the following relations 



d mi d d 



dxi M dX dx 
d m 2 d d 



(5) 



dx 2 M dX dx 
where M = mi + m 2 is implied. 

The eigenvalue equation, Hip = Wip may be written as 

jc \Jm\c 2 - h 2 Vi + c\jm\c 2 - k 2 V 2 2 - ip = (6) 

— * — * 

using pi = —ihVi and p 2 = —i%V 2 in Eqn.(l). Operating on the left by (H + W), we 
obtain, 

(W 1 + H 1 )ip = (7) 
with Wi = c 2 {m\c 2 - h 2 V 2 ) + c 2 (m\c 2 - h 2 V 2 2 ) - W 2 



and Hi = 2c 2 ^Jm 2 c 2 - h 2 V 2 ^m 2 c 2 - h 2 V 2 2 . 

To eliminate the square root operators in Hi, we operate on the left of Eqn.(7) by 
(Wi — Hi) and obtain 

(W 2 - Hi) V = 0. (8) 
The operator equation of the above is 



(ml-h 2 V*)\(?(ml-h 2 Vl) 2 -W^ -^(m^-^V^m^-^V 2 .) iP = 0. (9) 
The above equation simplifies to 

c 4 ( m 2 c 2 _ + c 4 ^2 _ ^ y 2^ + ^4 _ £ 4 (^2 _ ^2 y 2) ^^2 _ 

-2W 2 c 2 (m\<? - U 2 V 2 ) - 2W 2 c 2 (m 2 2 c 2 - h 2 vfj\ ip = 0. (10) 
Noting that Vi = -V 2 = V implies that V 2 = V 2 , = V 2 . 
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Ignoring the translational motion of the system and holding the centre of mass at 
rest, Eqn.(lO) reduces to 

W A +(m A c 8 +mtc 8 -2m 2 m 2 c 8 ) -2W 2 m\c A -2W 2 m\c A +4W 2 c 2 Ti 2 V 2 ] ip = 0. (11) 

It may well be checked that in the case of equal masses, mi = m 2 = m , the above 
equation reduces to 

which is in agreement with the work of Rao et.al 2 . 

The more general Eqn.(ll) may be written in an elegant form as 

y2 _|_ W 2 - 2 (m 2 c 4 + m^c 4 ) ^ (m\c A - m\c AV 



i> = (13) 



4c 2 k 2 4W 2 c 2 h 2 

To test the validity of this equation, we consider various cases. 

2.1 Single Particle Rotor 

For a single particle of rest mass mo, rotating about a fixed origin, Eqn.(6), with 
mi = m and m 2 = 0, would yield, 

^ + F^)* =o (l4) 

OR 

Expressing V 2 in spherical polar co-ordinates and noting that ~~q^ = 0, we write 

1 8 ( . n di>\ 1 d 2 i) (W 2 -e 2 \ , n 

7^ sin ^ + 1^—271 ir£ + — 553- W> = (15) 



r 2 sinfl 89 \ 89 j r 2 sin 2 fl 8(j) 2 \ c 2 h 2 
where e = m c 2 . 

With r = a, I = mod 2 , the moment of inertia of the particle about an axis passing 
through the origin, the above equation may be written as 

1 8 ( . n 84)\ 1 8 2 i) I(W 2 -e 2 ) , n 
smS <96> \ 96* / sin^ <90 2 eh 

Comparing this with the corresponding non-relativistic Schrodinger equation, we imme- 
diately obtain, 

^-(W 2 -e 2 )= 1(1 + 1), (17) 
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where / is an integer. 

It follows from above that the energy eigenvalues are 



W = ^ + '±^ (18) 

Non-relativistic limit : 

Taylor expansion of Eqn. (18) gives 

Defining W — e = E^ r , we get the well-known expression for the energies 

E Nr = /(/ Q + y 2 , where 1 = 0,1, 2, ... (20) 

2.2 Homonuclear system 

For a system of two particles of equal masses rotating about the centre of mass, 
77ii = m 2 — m o ) I — l^d 2 where the reduced mass [i = ^r and e = 2moc 2 , the rest 
energy of the system, Eqn. (13) transforms into 

„ 2 , (W 2 -e 2 )m . . 

V 2 ^ + 1 oft2 V> = (21) 

which is the appropriate Klein-Gordon equation of the system. As before the eigenen- 
ergies correspond to Eqn. (18). 

2.3 Heteronuclear System 

A more general case is a system of two non-interacting particles of unequal masses 
rotating about their centre of mass. The energy eigen-vlaues are obtained by solving 
Eqn. (13). 



Expressing V 2 in spherical polar coordinates, as before and observing that it trans- 
-L 



forms to the ( -f^- ) operator with eigenvalues 1(1 + 1), we write 



a 2 {W 2 - 2 (m 2 c 4 + m|c 4 )} a 2 (mfc 4 - m 2 c 4 ) 2 _ 

Writing ct = mic 2 , j3 = m 2 c 2 , it is obvious to check that the above equation may 
be written as 

a 2 W 4 — (A + B) W 2 + C = 0, (23) 
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where 



A = 2a 2 (a 2 + (3 2 ) 
B = 4l (l + l)c 2 h 2 
C = a 2 {a 2 - [3 2 f 

Eqn.(23) is solved using Mathematica 4 . Obviously, from the solution of the above 
equation for equal masses one should recover the energies of the homonuclear system. 
Hence of the four solutions, two of them being zero and the other being negative, are 
rejected and the only admissible solution is 



W= 



21 



(rr%c 4 +rr%c 4 )+ — (I +l)c 2 h 2 + — J(m 2 c 4 a 2 + c 2 h 2 l (I +l))(m|c 4 a 2 +1(1 +l)c 2 h 2 ) 



(24) 



The eigenenergies of a heteronuclear rotor are explicitly obtained from the above equa- 
tion 

Non-relativistic limit 

To arrive at the non-relativistic limit of Eqn.(24), we first square the expression and 
then use the Taylor expansion for the term within the square root. This leads to 



W 



, A 1(1 + 1) c 2 h 2 

1 + ^— n + 



(mic 2 - m 2 c 2 ) 2 1(1 + 1) c 2 h 
{mic 2 m 2 c 2 a 2 + 1(1 + 1) c 2 h 2 } 



.2*2 



(25) 



where e = mic 2 + m 2 c 2 is the total rest energy of the system. 



The relativistic correction is made evident by the binomial expansion of Eqn.(25) 
which gives 



2l(l + l)c 2 Ti 2 1 (m lC 2 -m 2 c 2 y 1(1 + 1) c 2 h 



2^ 



.2*2 



W 



e + 



a?e 



+ - 



2 e (mic 2 m 2 c 2 a 2 + 1(1 + 1) c 2 h 2 ^j 



(26) 



While the first term corresponds to the rest energy of the system, the second refers 
to the NR term and the third gives first relativistic correction. 



It is straightforward to check that 



E Nr — 



21(1 + l)h 2 n 
1 M 



(27) 



where /i = J^^^ ; I — ^ an d M = mi + m 2 as before. 



In the case of equal masses, we recover the well-known expression for the energies of 
the non-relativistic rigid rotator. 
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Relativistic Rotational Coefficient 

In spectroscopic calculations, the rotational constant, B, is an important parameter, 
characteristic of the molecule 5 . It is seen that the energy spectrum of a non-relativistic 
rigid rotator consists of a series of equidistant lines, the first one lying at 2B and the 
separation of successive lines also being 2B. 



Eqn.(26) may be written as 



W KG = e + l(l + l)hc 



h L (mi - m 2 ) 2 a 2 



An 2 Ic M I 4 ( mi m 2 c 2 a 2 + l(l + l)ti 



(28) 



Apparently, the second term within the brackets may rightly be called the 'energy 
correction term' resulting due to unequal masses. This term highlights the effect of going 
from homonuclear to heteronuclear species. 

The above equation may be put in a compact form as 

W KG = e+ 1(1 + 1) he B Rel (29) 



where 



_ h 2/i J K- rn 2 ) 2 a 2 

"** An 2 Ic M 1 4 ( mim2 c 2 a 2 + 1(1 + l)h* ' 



Bffei may be called the /-dependent Relativistic Rotational Coefficient. In the case of 
equal masses, 

Bb * = 8^Tc ' (31) 

which is in agreement with the well-known rotational constant, B of non-relativistic 
theory. Thus Eqn. (31) for homonuclear molecules is only a special case of the more 
general equation (Eqn. 30) applicable to heteronuclear molecules. 

Energy Spectrum 

Starting from Eqn. (25) and carrying out the Taylor expansion upto the first order, 
one can readily obtain 

W KG = e+ 1(1 + 1) hcB Rel - l H±}l^ki (32) 
which may be written as 

Wl= e+ l(l + l)hc(B + B l )- L ^ t + i; n ° (B + B{f (33) 

where 
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„ h 2u , , 

and 

5, = B , K^-^)^ 2 (35 ) 

4 (m lC 2 m 2 c 2 a 2 + Z (Z + 1) c 2 ft 2 ) 

Essentially, for a rotor having unequal masses, each level is shifted differently. The 
rotational constant is trivially modified and Bi may be interpreted as the relativistic 
correction term to the rotational constant, for a given transition. 

Eqn.(33) may also be written as 



;2 n , i \2 7 2 2 

Wi = Wq + 1(1 + 1) he Bi (Bf + 2BB t ) (36) 



) 2 h2„2 

2e 

where 



w^ t+ i ( i + i)hcB- (37) 

It is interesting to note that for equal masses, Bi = and W\ = Wo, the emperical 
formula obtained by Rao et.al. 

Wi+i-Wi 

Introducing the wave number, v\ = , we obtain for the transition I + 1 — > I, 

jy l = T 1 +T 2 + T 3 + T 4 + T 5 (38) 

where 

T x = 2 (Z + 1) B 

T 2 = J^B 2 (Z + l) 3 

T 3 = (/ + l){(/ + 2) B l+1 -lBi} 

T 4 = -^^{(l + 2fB? +1 -PB?} 

T 5 = -^^2B{(l + 2fB l+1 -PB l } 

A more rigorous calculation shows that the terms T 3 ,T 4 and T 5 are mass-dependent 
and are explicitly given by 

T _ B (/ + 1) m 1 c 2 m 2 c 2 (m x c 2 - m 2 c 2 f a 4 

3 ~ 2 (m lC 2 m 2 c 2 a 2 + (/ + !)(/ + 2) c 2 ft 2 ) (m lC 2 m 2 c 2 a 2 + /(/ + !) c 2 U 2 ) 



Ta 



B 2 hc (l + l) 2 (l + 2)m 1 c 2 m 2 c 2 (m 1 c 2 -m 2 c 2 ) 2 a 6 \m lC 2 m 2 c 2 a 2 -21 (Z + l)c 2 k 2 



lfc (m x c 2 m 2 c 2 a 2 + (Z +1) (Z +2) c 2 /* 2 ) (m lC 2 m 2 c 2 a 2 + 1(1+1) c 2 h 2 ) 



(40) 
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5 2 /ic (Z + 2) 2 (m lC 2 - m 2 c 2 ) 2 a 2 {2/ (/ + 1) c 2 ft 2 + (Z + 2) m lC 2 m 2 c 2 a 2 } 
2 e (m lC 2 m 2 c 2 a 2 + (Z + 1) (Z + 2) c 2 Zi 2 ) (m lC 2 m 2 c 2 a 2 + /(/ + !) c 2 k 2 ) 



(41) 



It is straightforward to check that T 3 ,T 4 and T 5 vanish when the masses are equal. 
While the first two terms of Eqn.(38) correspond to the usual expression for the wave 
number of the spectral line, the last three terms arising mainly due to unequal masses 
do contribute to the small relativistic correction in the wave number. 

Thus the application of relativistic quantum mechanics to the investigation of the 
spectra of heteronuclear linear diatomic molecules using the KG equation entails great 
mathematical complexities in all. 

First Rotational Line 

In the spirit of non-relativistic approximation, we retain T\ and T 3 and neglect all 
other terms in higher powers of B, to obtain 

v 2 = 2 (Z + 1) B + [l + 1) {(/ + 2) B l+1 - IBi} . (42) 



The first rotational line occurs at Z = 0, leading to 

u = 2B + 2B 1 . 



(43) 



It is trivial to check that 



B 
~2 



M 2 c 2 a 2 + 8h 2 



jiMc 2 a 2 + 2fi 2 _ 

where m\m 2 = /iM is implied. In units of Compton wavelength expressing a 



(44) 



and a 



h/jic 



, we obtain 



B 



v Q = 



a 2 + 8 
aa + 2 



h/Mc 



(45) 



For equal masses, &o — 2-B, by virtue of Eqn.(42). It is thus seen that the allowed 
transitions for homonuclear diatomic molecules are regularly spaced at an interval of 
2B. However, for a heteronuclear KG rotor, the lines are no longer exactly equidistant. 
Apparently, the measurement and identification of one of the spectral lines allows us to 
calculate the moment of inertia and then the bondlength of the molecule. 
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3 Results and Discussion 

A rigid rotator approximates a diatomic molecule when vibrational motion is ignored. 
In literature, the energy specturm of a rigid rotator is investigated non-relativistically 
using the Schrodinger equation. In the present article, we have introduced the notion of 
the 'relativistic rotator' and obtained the eigen energies relativistically. To investigate 
the energy spectrum of spin zero homonuclear and heteronuclear systems, we solve the 
Klein-Gordon equation in the centre of mass frame. Step by step elimination of the 
square root of an operator results in the fourth order equation in W, from which the 
eigen-energies are extracted. This novel approach can be extended to spin half systems, 
as well. 

Apparently, it is seen that the eigenenergies of heteronuclear systems reduce to those 
of homonuclear species when mi = m 2 , and have an approapriate non-relativistic limit. 
It is seen that the level spacing which is equidistant for a non-relativistic rotator is not so 
for a relativistic rotor. Interestingly, the first rotational line of a non-relativistic rotator 
and that of a Klein-Gordon rotor occurs at 2B. While it is known that the pure rotation 
spectrum of a diatomic molecule consists of equispaced lines, that of a relativistic rotator 
is not so. 

It has been established that the rotational constant of homonuclear and heteronu- 
clear diatomic molecules depend on the electronegativity of the constituent atoms 6 . 
Relativistic studies, apart from having immense pedagogical value are also of practical 
use in the light of ever increasing accuracies of modern experimental methods. 
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